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Abstract 



The morphological stability of two-dimensional islands nucleated on a sub- 



C^ I strate during vacuum or vapour-phase atom deposition is investigated. Using 

a. 

simple scaling arguments, it is shown that, contrary to expectation, dendritic 



islands may be converted into compact ones by increasing the deposition rate, 
provided that the size of the critical nucleus is large enough. Implications for 



1^ \ recent observations of Pt deposition on Pt(lll) are discussed. 
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Since the seminal work by Mullins and Sekerka |]l|, it is known that compact growth 
shapes are unstable, in general, when growth takes place by diffusion in some medium. 
In particular, scanning tunneling microscopy (STM) demonstrates the occurrence of fractal 
two-dimensional (2D) islands during metal-on-metal epitaxial deposition at low temperature, 
as reported by several authors |@-0. The physical mechanism leading to ramified — possibly 
fractal — growth shapes is the same as in the diffusion limited aggregation (DLA) model of 
Witten and Sander p. The deposited atoms (adatoms) diffuse on the surface, until they 
meet another adatom or an island. Any fluctuation of the shape leading to a protrusion 
tends to be amplified during growth, since diffusing adatoms find more easily the protrusion 
than the flat parts, which lag behind (c/. the effect of a point in electrostatics). At low 
enough temperature the interatomic bond cannot be broken, and island restructuring by 
detachment-reattachment of bound atoms is prevented. If diffusion of bound atoms along the 
island edge is also hindered, as it is the case for (111) metal faces, fractal shapes occur. This 
instability is intrinsic of diffusion limited growth at low temperature. Indeed, non-compact 
or, generally speaking, dendritic growth shapes appear at any temperature, provided that 
smoothing processes, e.g. edge diffusion, are not able to effectively counterbalance the 
intrinsic instability [^. 

On the theoretical side, a detailed analysis of the morphological stability of a 2D island 
has been performed by Avignon and Chakraverty [0 for small adatom supersaturations. The 
regime of high supersaturation, when the size of a critical nucleus is 1 [^ and classical fractals 
appear, has been investigated by Pimpinelli et al. 0, with an emphasis on the stabilizing 
action of several matter-transport processes (edge diffusion, detachment-attachment). The 
same situation has been subsequently addressed with Kinetic Monte Carlo simulations by 
Bartelt and Evans [|1^], who studied the shape transitions of a single island, and by Bales and 



Chrzan |^, who studied the compact-fractal transition during epitaxial deposition. Both 
works assume edge diffusion as the smoothing mechanism, and forbid atom detachment 
from the island. Their results are in good agreement with the analytic computations in 0. 
Also, the dependence of the fractal instability on the geometry of the substrate has been 



investigated by Zhang et al. [jT2|. 



One of the resuhs of such studies is that, not surprisingly, the fractal instability is 
favoured by a high deposition rate: if the time interval between the arrival of two adatoms 
onto the island edge is short, edge diffusion cannot exert its smoothing action, because before 
the first atom has had the time to search for a site of high coordination — which is a condition 
for a compact shape — it is blocked by the second incoming atom — and limited edge mobility 
clearly favours the instability. It appears thus obvious that such morphological instabilities 
are due to the far-from-equilibrium nature of the deposition process, and that going farther 
out of equilibrium (growing faster) can only make things worst. 

However, during atom deposition on a substrate, many islands are simultaneously nucle- 
ated, which then grow laterally and compete for capturing additional incoming atoms. As 
a result, each island has a "capture area" limited by its neighbours and its advancing edge 
slows down in approaching the capture area of other islands. A sort of repulsion between 
island edges appears, which acts as an indirect shape-smoothing mechanism. In general, 
increasing the deposition rate increases the island density, and thus the effective repulsion. 
Indeed, I will show that, under appropriate conditions, increasing the deposition rate makes 
the islands more compact. 

Let me assess the picture by recalling known results. Note that throughout this Letter, 
lenghts are measured in units of the interatomic distance a, so that a = 1. 

The evolution of the island morphology is determined by the balance between destabiliz- 
ing and smoothing processes. Since kinetics is involved, it is a matter of timescales. Calling 
F the deposition rate, and i"^ the average size of the capture area of an island, the average 
number of atoms reaching the island per unit time is $ = Fi"^. The time elapsed between 
two atoms successively hitting the same site on the island edge is then r = 1/$. For the 
instability to be checked, the time r must be larger than the average time taken by an atom 
to diffuse along the island edge. If Dg is the edge diffusion constant, and L^ the average 
island size, there is no instability as long as r > L'^/Df,. Thus, the instability takes place 
when L^ exceeds the critical value L^ given by |]9|JT0[ 



Ll ^ DJ^ . (1) 

Therefore, the instabihty has in this case a threshold which decreases with increasing 
F. Note at this point that confusion should be avoided between the critical size for the 
appearance of the instability, L^, and the size of the critical nucleus, n. The latter is the size 
beyond which an island grows continuously without net loss of atoms. Aggregates smaller 
than the size n may also shrink to smaller sizes (see note p|). The DLA-type situation 
corresponds to n = 1 and very low temperature. At higher temperature, the possibility of 
atom emission from kinks has to be accounted for, and the edge diffusion constant De must 
be replaced with a kinetic coefficient Fg, which also depends on the formation energy of a 
kink, and thus on the step line tension (see Refs. |l|,§] and point (iii) below.) 

When just one island is present, P' is just the substrate area, and $ is just the number 
of deposited atoms per unit time. Equation (|l]) then states that Li varies as 1/-F. This is 
an extreme case. Indeed, when many islands are nucleated on the surface, as it is usually 
the case in epitaxial deposition, the capture area f' has to be computed self-consistently 
by accounting for the competing presence of the other islands. This task needs considering 



the details of island nucleation. It has been performed in a number of papers |T3|-p3|, and 
the result is the following: C."^ is proportional to a power of the ratio Ds/F, where Dg is the 
surface adatom diffusion constant, the power being a function of the critical nucleus size n: 

f ^ (D,/F)"/(2+-) . (2) 

Since each island sits, on average, at the center of a region whose area is £^, the average 
island density is A^ = l/£^. Hence, the average island size L^ and the total deposited 
matter — the surface coverage 6 = Ft — are related via L'^N = 6, or 

L^ = 9f (3) 

Using $ = Fi"^, and equations (|1|) and (|^) the critical island size L^ can be converted 
into a critical coverage 6c as follows 



Oc ^ Ll/e ^ (TJD,) -^ . (4) 

Replacing £^ with its expression (g) yields 

e, ^ ivjD,) (F/D,)("-2)/("+^) . (5) 

This result is new. 

When n = 1, equation (^ becomes 

d, ^ {TJD,) {DjFf . (6) 

Equation (^) agrees well with Bales and Chrzan's simulations [|rT|. It shows that for n = 1, 
the critical coverage for the appearance of the instability decreases with increasing F as 
expected. However, we immediately see that, contrary to expectation, when n > 2 the 
opposite is true: the critical coverage in (j^) increases with increasing F. The increased 
growth rate stabilizes the compact shape! 

Other smoothing mechanisms can be investigated by this approach. One only needs to 
know for each healing process, the typical relaxation time of a perturbation of wavevector 
q. These timescales are contained in Ref. |T^. Instead, I rely here on an argument a la 
Mullins and Sekerka |1|,|3J^; which makes use of more macroscopic concepts as the curvature 
of the island profile, and the excess chemical potential of the curved parts (Gibbs-Thomson 
relation [|T3). Indeed, consider a circular island of radius R. Let 5q be the amphtude of 
a deformation of the circular shape whose wavevector is q. Then, its rate of variation, 5q, 
results from the destabilizing diffusion contribution [|I], 



D,Vc ■ n ^ v,tcpq5g ~ qFf/R Sg , (7) 

where c is the adatom density, n the local step edge normal and fstep ~ F(? / R is the ad- 
vance velocity of the island edge step, plus the various smoothing terms. Three of them will 
be considered here: i) exchange of atoms between the island and a 2D adatom gas on the 
substrate; ii) atom detachment, diffusion on the surface and reattachment; iii) edge diffu- 
sion again, for completeness. Smoothing occurs because, according to the Gibbs-Thomson 



relation [|T5|, an increase in the curvature of an interface raises the chemical potential of the 
curved part, by a quantity 6fi ~ 'yq'^Sg, where 7 is the interface stiffness, and q the wavevec- 
tor of the deformation. The density of atoms near the interface is accordingly increased with 
respect to the equilibrium density Ceq. 

Cexcess = Cgq eXp(/3 S fx) ^ Ceq(l + P^Q^Sg) , (8) 

where l/P = ksT. 

i) When the excess atoms detach from the island and exchange with a 2D adatom gas 



on the substrate at rate A;exch, the smoothing term has the form |0 ^exchCexcess ~ rexchQ'^'^q, 



where (|) has been used, and Fexch = ^exchCeq/97- It is then straightforward to show in this 
case that no critical island size must be exceeded for the instability to appear. Any size is 
unstable, provided that F/Dg > (Dg/TexY'^'^'^^^'^- If the opposite is true, no size is unstable. 
When this smoothing process dominates, a "typical" situation occurs, in which increasing 
F is destabilizing. What smoothing process dominates is a system-dependent matter, of 
course. 

ii) When the excess atoms detach from the island and surface diffuse towards places of 
lesser curvature, the smoothing term has the form 

D^VCexcess ■ n ^ T sQ X q^Sq (9) 

where (H) has been used, and Tg = DgCeqP'j. Equating (|^) to (0) the instability threshold 
Qc is found, 

Rq, ^ {FfR/r,y/' . (10) 

The wavevelength 1/q of a shape perturbation cannot be arbitrarily large. Indeed, if R is 
the island radius one has 1/q < R. The critical radius Re can then be found by letting 
q = l/R, which yields Re ~ Ts/{Fi'^), and eventually leads to 

Oe ^ iVjD^f (F/Z},)("-^)/("+2) . (11) 



when (H) and (H) are used. The quahtative result is the same as that of (j^), i.e. when n is 
large enough (ra > 4 here), increasing F increases the stability of the compact shape. 

iii) When excess atoms diffuse along the interface, they give rise to a smoothing term of 
the form ^ 

DeV^Cexcess ~ ^ eQ^ X q'^Sg (12) 

where (|) has been used, and Fg = DeCeqP^. Equating ([T2|) to (0) the instability threshold 
Qc is found, 

Rq,^{Fe^RyT,y/\ (13) 

The critical radius Re can then be found by letting q = 1/R, which yields (|l|), and eventually 
leads to (^ when (^ and (|^) are used. The reader can check that the alternative method 
based on Ref. [|1^ gives the same results in all cases. 



The interesting point is that the scenario sketched here appears to have been observed. 
Michely, Hohage and Comsa [^ report that the triangular islands formed during Pt deposi- 
tion on Pt(lll) at 445 K and F = 2.1 x 10~^ ML/s, are morphologically unstable. Under 
these conditions, they exhibit dendritic shapes at a coverage 6 = 0.1 monolayer (ML). On 
increasing the flux by two orders of magnitude and leaving all other conditions unchanged, 
the islands appear to be quite compact. Stability is indeed obtained by depositing faster! 
A quantitative comparison of the qualitative arguments given in the present Letter and ex- 
periment is obviously out of the question. Pt(lll) is known as a system with complicated 
behaviour PJ^, which only detailed kinetic Monte Carlo simulations can hope to grasp 
|l^,|T8|. I could only venture to suggest that for such a system at a relatively high substrate 



temperature, the critical nucleus size is likely to be large, possibly n = 6 |[T^, so that both 
mechanisms (ii) and (iii) would predict the observed stabilization at increasing deposition 
rate, if they happen to dictate step smoothing on Pt(lll). 
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